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In this paper we introduce a hands-on activity in which introductory astronomy students act as 
gravitational wave astronomers by extracting information from simulated gravitational wave signals. 
The process mimics the way true gravitational wave analysis will be handled by using plots of a 
pure gravitational wave signal. The students directly measure the properties of the simulated signal, 
and use these measurements to evaluate standard formulae for astrophysical source parameters. An 
exercise based on the discussion in this paper has been written and made publicly available online 
for use in introductory laboratory courses. 



I. INTRODUCTION 

Observational astronomy stands at the threshold of an 
era where gravitational wave detectors are a tool which 
regularly contributes important information to the grow- 
ing body of astrophysical knowledge^ Ground based de- 
tectors such as the Laser Interferometer Gravitational- 
wave Observatory^ (LIGO) and the forthcoming space 
based detector the Laser Interferometer Space Antenna^ 
(LISA) will probe different regimes of the gravitational 
wave spectrum and observe sources that radiate at differ- 
ent gravitational wavelengths. Unlike their cousins, tra- 
ditional electromagnetic telescopes, gravitational wave 
detectors are not imaging instruments. How then does 
a gravitational wave astronomer take the output from a 
detector and extract astrophysical information about the 
emitting sources? This paper introduces a hands-on ac- 
tivity in which introductory astronomy students answer 
this question. 

Traditional astronomy is often presented through the 
medium of colorful images taken with large scale tele- 
scopes. In addition to studying images, astronomers 
learn about astrophysical systems by collecting data at 
multiple wavelengths, using narrow band spectra, mea- 
suring time varying light curves, and so on. It is often 
the case that the core physics governing the evolution of 
these distant systems is deduced from the physical char- 
acter of the observed electromagnetic radiation, rather 
than from the imagery that is used to illustrate the sci- 
ence for other audiences. 

Gravitational wave astronomy is completely analogous 
to its electromagnetic cousin, with one important distinc- 
tion: there will be no image data, because gravitational 
wave detectors are not imaging instruments. Gravita- 



tional wave observatories like LIGO and LISA return 
a noisy time series that has encoded within it gravita- 
tional wave signals from one or possibly many overlap- 
ping sources. To gain information about the systems 
emitting these gravitational wave signals requires the use 
of time series analysis techniques such as Fourier trans- 
forms, Fisher information matrices, and matched filter- 
ing. Recently an activity has shown how students can 
emulate the match filtering process by comparing ideal 
signals to mocked noisy detector output^ In this paper, 
a procedure is described whereby students can analyze a 
simulated gravitational wave signal and extract the astro- 
physical parameters which describe the radiating system. 
The goal is to introduce students to how gravitational 
wave astronomers learn about sources of gravitational ra- 
diation in a fashion suitable for classroom or laboratory 
exercises related to this modern and emerging branch of 
observational astrophysics. 

The rest of this article introduces some basic back- 
ground of gravitational radiation and then describes the 
activity. Section [Q] outlines the theory connecting the 
structure of gravitational waves to astrophysical param- 
eters, and Section IITT1 illustrates the characteristic wave- 
forms from a typical binary system. Section llVl illustrates 
a procedure where measurements made from waveform 
plots, together with the theory of waveform generation, 
can be used to extract the astrophysical parameters (or- 
bital period, distance, etc.) of the system emitting grav- 
itational radiation. Section Ivl discusses implementations 
and extensions for this activity in an introductory as- 
tronomy course. The analysis described in the paper has 
been implemented in an activity format, complete with a 
keyed solution for the instructor, and is publicly available 
online. 5 
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II. GRAVITATIONAL WAVE PRODUCTION IN 
BINARIES 

In electromagnetism radiation is produced by an ac- 
celerating charged particle. Similarly in general relativ- 
ity gravitational radiation is produced by an accelerating 
mass. To be precise, gravitational waves are produced by 
a time varying mass quadrupole moment. The reason for 
this is straightforward^ Monopolc radiation is prevented 
due to conservation of mass, while dipole radiation does 
not occur due to the conservation of momentum. This 
leaves the quadrupole as the leading order term in the 
multipolc expansion of the radiation field. A simple and 
common example of an astrophysical system with a time 
varying quadrupole moment is a binary star system. 

For a circular binary system, where the components 
are treated as point-like particles, the gravitational wave- 
forms take on the seductively simple form 



h{t) = .4(f)cos($(f)) , 



(1) 



where h(t) is the gravitational waveform (also referred 
to as the gravitational wave strain), .4(f) is the time de- 
pendent amplitude, and <f>(f) is the gravitational wave 
phase. The amplitude A(t ) can be expressed in terms of 
the physical parameters characterizing the system, 

where G is Newton's gravitational constant, c is the speed 
of light, r is the luminosity distance to the binary, and 
Pgw{t) is the gravitational wave period. The quantity 
M = (MiM 2 ) 3/s (Mi + M 2 )~ 1/5 is called the chirp mass 
and appears repeatedly in gravitational wave physics, 
making it a natural mass scale to work with. The origin 
for this nomenclature will become evident shortly. The 
waveform phase $(f) is given by a the integral: 



$(f) = $ + 2tt / 
Jo 



d£ 



Pgw(t') 



(3) 



where $o is the initial phase. 

As in electromagnetism, gravitational waves have two 
independent polarization states. For a binary system the 
two states are related by a 90° phase shift. Consequently, 
Eq. ||TJ captures the functional form for both polariza- 
tion states. For the purposes of this paper only a single 
polarization state and its associated waveform will be 
discussed. 

Gravitational waves carry energy and angular momen- 
tum away from the binary system causing the orbital 
period to decrease with time according to± 

P orb (t)= (p 8/3 -^ 3/8 , (4) 

where P is the orbital period at time f = 0, and k is an 
evolution constant given by 
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As a consequence of the ever shortening orbital period, 
the two binary components will slowly inspiral, eventu- 
ally colliding and coalescing into a single remnant. Under 
the assumption of point-like particles made here, this for- 
mally occurs when P or b{t) = 0. 

Note that Eq. @ gives the orbital period not the 
gravitational wave period. Careful scrutiny of Eqs. (J5J) 
and will reveal that the gravitational wave period 
Pgw(t) is the quantity which appears in the description 
of the waveform. Fortunately, for circularized binary sys- 
tems, P rb(t) and P gw (t) are simply related: 



Porb(ty ^Pgwify 



(6) 



(5) 



The simple factor of two stems from the fact that the low- 
est possible order for gravitational radiation production 
is the quadrupole order. Moreover, quadrupole moments 
are invariant under a 180° rotation, yielding a factor of 
two per complete orbit. 



III. WAVEFORMS FROM A BINARY SYSTEM 

As an illustrative example of the kind of waveforms we 
expect from binaries, consider a binary neutron star sys- 
tem with Mi = M 2 = 1.4 Mq (M = 1.22M ) located 
at the center of the galaxy r = 8 kpc away. For the 
activity we will consider waveforms generated at two dis- 
tinct times in the binary's evolution. The first waveform 
we will consider is ~ 10 6 years before coalescence. Dur- 
ing this phase the gravitational wave frequency is in the 
regime that will be detectable by the spaceborne LISA 
observatory, which has a principle sensitivity in the range 
of 1CP 5 Hz to 1 Hz. The second waveform considered will 
be during the final second before the neutron star binary 
coalesces. The gravitational wave frequencies during this 
phase are in the regime that will be detectable by the 
terrestrial LIGO observatory, which is sensitive to gravi- 
tational wave frequencies between 10 Hz and 10 3 Hz. 



A. Far from Coalescence 

Figure^shows the emitted gravitational radiation long 
before the binary components coalesce. During this era 
of the binary evolution, the gravitational waves are es- 
sentially monochromatic; the orbital period is evolving 
too slowly to detect a frequency derivative term. 

For monochromatic signals like this, the only measur- 
able properties of the gravitational waveform are the pe- 
riod P gw (and the orbital period P or b through Eq. JBJ), 
amplitude A, and initial phase 4>o. Even though the 
waveform equations depend on the chirp mass M. and 
the luminosity distance r, it is not possible to solve for 
their values from the data provided by the monochro- 
matic waveform. Not enough information exists to com- 
pletely solve Eqs. J5J) and (0} together for both quanti- 
ties. This can be seen by considering the relative size 
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FIG. 1: The gravitational waveform for a binary system con- 
sisting of two neutron stars far from coalescence and located 
at the center of our galaxy. 

of the two terms in Eq. using the binary neutron 
star chirp mass A4 = 1.22 M Q , it should be evident that 
the second term is completely negligible compared to the 
period Pn of the wave shown in Fig. ^ In the parlance 
of gravitational wave astronomy, there is a mass- distance 
degeneracy in the waveform description, analogous to the 
familiar mass-inclination degeneracy in the electromag- 
netic observations of spectroscopic binaries. This degen- 
eracy is a well known problem, but as the next section 
shows, it can be broken if the orbital period of the binary 
evolves during the gravitational wave observations. 



B. Near Coalescence 

Inspection of Eq. (0J shows that as time goes on, the 
emission of gravitational waves causes the orbital period 
to grow shorter, and as a result the frequency of the emit- 
ted waves increases. Similarly, consideration of Eq. @ 
shows that as the wave period decreases, the time de- 
pendent amplitude A(t) increases. This is characteristic 
behavior for gravitational waves emitted just prior to a 
source coalescence, and is known as a chirp. The chirp 
waveform emitted by the example binary neutron star 
system just prior to coalescence is illustrated in Fig. [21 

Any binary signal which evolves appreciably during the 
gravitational wave observation is called a chirping binary. 
In these cases, the mass parameter M. which appears in 
the amplitude A{t) and in the period evolution constant 
k can be determined from measurements of the evolving 
signal. For this reason, the mass M. is called the chirp 
mass. To leading order in gravitational wave production, 
it is not possible to measure the individual masses, only 
the chirp mass. Consequently, it is not possible to dis- 
tinguish between binaries with the same chirp mass. For 
example, the binary neutron star considered in this pa- 
per with Mi = M2 = 1.4 M Q has roughly the same chirp 



- 






































J 


lllifflll 

1 ( 1 




M| | 


Piiiim 


■Ill 


HI 


i A? 

in 


w 








lllllllllll 


ii'liiili 


if 


Si, 1 , 1 ,!,'] 


III 


•- 

umi 


I ill 1 !! 

II i' "'i 


TT 

S 
'ii ||i| 


|1 
"III 




1 









































0.2 0.4 0.6 0.8 1 

t (seconds) 



FIG. 2: The waveform over the last second before coalescence. 
Since the signal's amplitude and frequency is increasing with 
time, these types of systems are said to be chirping. 



mass as a binary with an Mi = 10 M Q black hole and a 
M 2 = 0.3 M Q white dwarf. 

To extract the chirp mass from measurements of the 
gravitational waveform, consider two small stretches of 
the chirping waveform. Figure shows the waveform 
from s < t < 0.05 s, and Fig. 0] shows the waveform 
from 0.9 s < t < 0.92 s. The waveform is appreciably 
different between these two snapshots, both in amplitude 
A{t) and in period P gw (t). This allows the degeneracy 
found in the monochromatic signal case to be broken, 
because the gravitational wave period can be measured 
at two different times and used in Eq. (0} to solve for the 
chirp mass M.. 



IV. MEASURING GRAVITATIONAL 
WAVEFORMS 

This section illustrates a procedure at the introductory 
astronomy level where students can make direct measure- 
ments from the figures in Section lTTTl using a straight edge 
and the axis labels. Using their measured data together 
with the theory presented in Section [HJ the astrophys- 
ical character of the system emitting the gravitational 
waveforms can be deduced. 



A. Monochromatic Waveforms 

Limited astrophysical information can be extracted di- 
rectly from Fig.Q] as will be the case with true monochro- 
matic signals detected by gravitational wave observato- 
ries. With limited assumptions more detailed informa- 
tion can be deduced, which will be valid so long as the 
assumptions are valid. A suitable extraction and analysis 
procedure for an introductory astronomy student would 
proceed in the following manner: 
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• The gravitational wave period P gw can be mea- 
sured directly from the figure. Since the signal is 
monochromatic, the binary is circular and the or- 
bital period P or b is obtained directly from P gw us- 
ing Eq. iJSJ. For the waveform in Fig.^cajcful mea- 
surement should yield a value of P gw = 1000 sec. 

• The amplitude A and the initial phase $o can also 
be measured directly from the figure. As noted 
in Section IHIAI no astrophysical information can 
be extracted from the amplitude alone. The initial 
phase is a simple quantity to measure, but does not 
represent any intrinsic property of the binary; its 
value is solely a consequence of when the gravita- 
tional wave observations began. To illustrate this, 
imagine relabeling the time axis in Fig. ^ to rep- 
resent a new observation which started somewhat 
later than the observation shown. The initial phase 
will have some new value, but the waveform itself 
does not change because the intrinsic properties of 
the binary did not change. 

• If a gravitational wave astronomer were to assume 
that the binary was a pair of neutron stars, the 
component mass values could be assigned as part 
of the assumption. Most neutron star masses clus- 
ter around M = 1.4 M Q , so a good base assump- 
tion is that each component of the binary has this 
mass. As noted in Section IIII Al this assumption 
can be a dangerous one, since similar chirp masses 
M. can result from significantly different systems. 
Other information, not present in the gravitational 
waveform, may help an astronomer feel more con- 
fident about such an assumption. For example, an 
associated simultaneous electromagnetic signal or 
the location of the source on the sky may favor one 
model of the binary over another. 

• If the masses are assumed, the orbital separation 
of the binary components, R, can be computed 
from the measured orbital period by using Kepler's 
Third Law: 

G(mi + m 2 ) = i? 3 . (7) 

\ + orb J 

For this example, the orbital separation is R = 
1.4 x 10" 3 AU = 2.1 x 10 8 m, or a little less than 
the separation of the Earth and the Moon. 

• If the masses are assumed, the distance to the bi- 
nary can be computed from Eq. J5J) and the mea- 
sured amplitude. If careful measurements have 
been made, the answer should be close to the value 
r = 8 kpc = 2.5 x 10 20 m. 

• Lastly, if the masses are assumed, it can be quanti- 
tatively shown that the monochromatic descriptor 
is a good one for this wave by computing the value 
of the second term in Eq. (0} and showing that it 
is negligible compared to the measured period Pq. 
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FIG. 3: The chirping waveform one second before coalescence. 
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FIG. 4: The chirping waveform one-tenth of a second before 
coalescence. 

B. Chirping Waveforms 

In the case of a chirping waveform, additional as- 
trophysical information associated with the system can 
be extracted directly from measurements of the wave- 
form without making underlying assumptions like those 
needed when the system was far from coalescence. To ex- 
tract information from the chirping waveform shown in 
Fig. El the two zoom-ins of the waveform shown in Figs. [31 
and0]will be used. A typical extraction procedure might 
look like this: 

• For each of the Figs. |3 and 0] measure the period 
of one cycle of the wave, and note the time t at 
which the periods were measured. The amplitudes 
A(t) should be measured for the same cycle as the 
periods. 

• If the period measured at time t± in Fig. [3] is Pq, 
and the period measured at time t 2 in Fig. 0] is 
Pgw (t) at time t = t 2 — 1\ , then Eq. Q can be used 
to deduce the chirp mass M. of the system. 
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• Once the chirp mass M has been determined, the 
distance to the binary can be computed by using 
Eq. (J2J) with the measured amplitude A(t) and pe- 
riod P gw (t) of each waveform. The results from the 
two figures can be averaged together to obtain a 
final result. 



V. DISCUSSION 

This paper introduced the core calculations in 
gravitational wave astrophysics an introductory as- 
tronomy student can perform in a laboratory set- 
ting to glean information about an astrophysical 
system. To compliment this article we have also de- 
veloped a student activity sheet and corresponding 
teacher's guide related to the exercises described in 
sections II V Al and II V Bl The complimentary ma- 
terial is available, along with a template activity^ at 
http : //cgwp . gravity . psu . edu/ outreach/ activities/ 

The activity described here is a simple introduction to 
how gravitational wave astronomers extract astrophysi- 
cal information from observed binary waveforms. Real 
signal analysis is a more complex endeavor than what 
has been presented here. The most significant challenge 



in the case of true data is identifying the signal buried in 
a noisy data stream. A common approach to this prob- 
lem in gravitational wave astronomy is to use template 
matching, which has been explored in a separate activ- 
ity^ If a signal is present in a noisy data stream, the 
template provides a way to subtract the noise away and 
leave a clean waveform behind. This is the assumed start- 
ing point in the activity developed here. Its from clean 
waveforms that astronomers will estimate the values of 
astrophysical parameters describing a source of gravita- 
tional waves. By sequencing the two activities a student 
is exposed, at least in an idealized way, to the methods 
used by gravitational wave astronomers to extract astro- 
physical information about the emitting systems. 
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